VORTICES AND SPONTANEOUS SYMMETRY BREAKING 
IN ROTATING BOSE GASES 
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Abstract. We present a rigorous proof of the appearance of quantized 
vortices in dilute trapped Bose gases with repulsive two-body interac- 
tions subject to rotation, which was obtained recently in joint work with 
Elliott Lieb [1]. Starting from the many-body Schrodinger equation, we 
show that the ground state of such gases is, in a suitable limit, well 
described by the nonlinear Gross-Pitaevskii equation. In the case of ax- 
ially symmetric traps, our results show that the appearance of quantized 
vortices causes spontaneous symmetry breaking in the ground state. 



1. Introduction 

In recent remarkable experiments [21 El HI IS] , the appearance of quantized 
vortices in the ground state (and low temperature equilibrium states) of 
rotating dilute Bose gases was beautifully demonstrated. These quantized 
vortices are a consequence of the superfluid nature of the system under 
investigation. In particular, since the system is almost completely Bose 
condensed, it behaves like a single quantum particle. 

The state of ultracold dilute Bose gases is usually described by means 
of the Gross-Pitaevskii (CP) equation [3 Ei El El lO]. This non-linear 
Schrodinger equation originates as the variational equation from the cor- 
responding GP energy functional, given by 

= mo\<l^)+^^9 [ \H^td\ . (1) 
Jm.3 

Here, cj) G L^(M^), and Hq denotes the one-particle Hamiltonian, describing 
the kinetic, potential and rotational energy of the particles. In fact, if Q 
denotes the angular velocity vector and V{x) the trap potential, Hq is, in 
appropriate units, given by 

Ho = -A + V{x) -n-L, (2) 

where L = — ixA V denotes the angular momentum operator. The parameter 
5 in ([T|) is nonnegative and measures the interaction strength among the 
particles. The trap potential V{x) is assumed to be locally bounded and to 
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increase fast enough at infinity in order to have the particles confined to the 
trap (and, in particular, to ensure that Hq is bounded from below). More 
precisely, we assume that 

lim (v{x)-hnAx\A =+oo. (3) 

|a;|— *oo \ A J 

Since —A — Q ■ L = (— -iV + A — |il A this implies the desired 

property. 

The GP energy is the minimal value of £''■'^[0] among all (appropriate 
normalized) functions (j), i.e., 

i?GP(5,f^)= inf ^GP^. 

\m\2=i- 

Using ([3]) and the fact that g > 0, it is in fact not difficult to show that the 
infimum is actually a minimum (see [H]). That is, there exists a minimizer 
of the GP functional ([T]). Note that, in general, there may be many different 
minimizers. In any case, any minimizer satisfies the GP equation 



- A<t){x) + V{x)4>{x) - n ■ L<p{x) + 8TTg\<f>{x)\^(t>{x) = fi(t>{x) 



where fi = E'^^{g,il.) + Ang \(j){x)\'^d^x is the corresponding chemical 
potential. 

For axially symmetric V{x), i.e., in case V{x) commutes with ■ L, the 
GP functional is invariant under rotation about the axis. It turns out 
that for any 0,^0, this rotational symmetry is broken in the GP minimizer 
for large enough interaction strength g 1 12\ I13j . This symmetry breaking 
is the result of the appearance of quantized vortices since, in case of more 
than one vortex, they cannot be arrange in a symmetric way. Note that, in 
particular, this implies that there will be many GP minimizers (for g large 
enough) . 

We remark that the phenomenon just described is a special feature of 
rotating systems and cannot be observed in a non-rotating system. In fact, 
for $7 = there is always a unique minimizer of the GP functional [llj . 

It turns out that the appearance of quantized vortices, and the resulting 
symmetry breaking, which we have just described, are not merely a property 
of the GP theory, but can actually be derived out of the underlying (many- 
particle) Schrodinger equation. This was proved in [1]. In the following 
sections, we will give a summary of these results, and we will explain the 
key ideas leading to their proof. 

2. The Schrodinger Equation for Many Particles 

Consider a quantum-mechanical system of a large number, A^, of bosons, 
with one-particle energies described by Hq (given in ([2]) above). We assume 
that the particles interact via a repulsive pair interaction potential Va{x). 
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The Hamiltonian for this system is given by 

N 

Hn = Y,h!,'^+ Y1 Va{x^-Xj), (4) 
i=l ^<i<j<N 

where the superscript {i) refers to the fact that Hq acts on the i'th variable. 
Since the particles under consideration are bosons, the Hamiltonian Hi\i 
acts on the subspace of totally symmetric functions in (^^L^(R'^), which 
we denote by Hn- 

The interaction potential Va{x) is assumed to be nonnegative and of short 
range. More precisely, it is assumed to have finite scattering length [111 I14j . 
denoted by a, which means that it has to be integrable at infinity (i.e., it 
has to decay faster than |x|~^). A typical example would be a hard sphere 
interaction, which formally means that Va{x) = oo for |x| < a and Va{x) = 
otherwise. We shall, in fact, choose some fixed (nonnegative) interaction 
potential w{x) with scattering length 1 and obtain Va{x) by scaling as 

Va{x) = a~'^w{x/a) . 

It is then easy to see that Va{x) has scattering length a. Moreover, a now 
appears as a parameter in the Hamiltonian if at, which can be freely varied. 
In particular, we can (and will) let a depend on A^. We note that this scaling 
of Va{x) is, of course, mathematically and physically equivalent to scaling 
the trap potential V{x) (and the angular velocity O) in an appropriate way, 
while keeping the interaction potential fixed. 

2.1. Ground State Energy. For fixed w{x) and V{x), we shall denote the 
ground state energy of Hn as Eo{N,a,ft), i.e., 

E„(N.a.a)= i„f 

Since the ground state energy per unit volume of a homogeneous Bose gas 
with interaction Va{x) at density p is given by Anap^ for low density [15j, it 
is reasonable to expect that Eq(N, a,Q) ~ NE'^^ {Na,Q) for dilute gases. 
Here, dilute means that a^p <C 1, where p denotes the average density. This 
condition is, in particular, satisfied if A ^ 1 and Na = 0{1). We call this 
the GP limit. In this limit, we have the following result [1]. 

THEOREM 1. For any g > and QeR^, 

(5) 



N^oo N ^ ' 



That is, for large A and a = 0(1/A), the ground state energy per particle 
is given by the GP energy with coupling parameter g = Na. Theorem [T] 
holds for all angular velocities 0, (satisfying the stability criterion ([3])). It 
extends previous results in the nonrotating case = [11]. 

Note that the right side of ([5]) is independent of the choice of the unsealed 
interaction potential w{x). In the dilute limit considered here, only the 
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scattering length a matters, and not the details of the interaction potential. 
Note also that the result cannot be obtained by simple perturbation theory; 
in fact, the / |i?i>|^ term in the GP functional is partly kinetic energy, and 
not the average value of Va{x) (which might even be zero, as in the case of 
the hard-sphere interaction). 

As will be pointed out in Subsect. 12.31 below, it is essential to restrict to 
symmetric wave functions (bosons) in Theorem [T] For the absolute ground 
state energy (defined as the infimum of over all wavefunctions, not 
necessarily symmetric ones), the result is wrong, in general. For the absolute 
ground state, the right side has to be replaced by minimizing a density- 
matrix functional instead |13j . 

2.2. Bose-Einstein Condensation. The GP energy functional ([I]) and its 

minimizers contain information not only about the ground state energy of 
the many-body Hamiltonian but also about the ground state or, more 
precisely, its reduced density matrices. Recall that for any wavefunction 
^ G T~Cn, its reduced one-particle density matrix is given by the kernel 

J^j^Hx, x') = N / ^{x,X2, ■ ■ ■ ,Xn)^*{x' ,X2, . . . , XN)d\2 ' ' ' d^N ■ 

Note that this defines a positive trace class operator on the one-particle 
space L2(]k3). 

The one-particle density matrix of a state ^ contains all the information 
about the system concerning expectation values of one-particle operators. 
It particular, the concept of Bose-Einstein condensation (BEC) is defined 

in terms of 7^^. 

Note that if ^' is normalized, i.e., ||^||2 = 1, then the trace of 7^^ is 

N. BEC means that 7^^ has an eigenvalue of order N. The corresponding 
eigenfunction is called the condensate wave function. For dilute systems, as 
we consider here, one expects in fact complete BEC, meaning that 7^^ is 
approximately a rank one projection, or 'y^\x,x') ~ N(j){x)(j){x') for some 
normalized (p G L^(M'^). 

In the non-rotating case = 0, complete BEC in the ground state of 
H]\f was proved in [16J. Moreover, it was shown that the condensate wave 
function equals the GP minimizer. Recall that in the case $7 = there is a 
unique minimizer of the GP functional ([T]) (up to constant phase factor, of 

course), which we denote by 0*^^. That is, if 7^^ denotes the one-particle 
density matrix of the ground state ^ of for = 0, then 

(6) 



hm l^«(:c,xO=0^P(x)</>GP(x') 



in the GP limit ^ 00, Na g. To be precise, the limit ([6]) holds in trace 
norm sense. Note that although a is scaled to zero in the limit considered, 
the right side of ([6]) depends on g = Na via 0*^^. 
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The corresponding result for Q ^ is necessarily more complicated be- 
cause of non-uniqueness of the GP minimizer (f)'^^ . It is actually more nat- 
ural to not just look at a ground state of Hn (which may not be unique in 
the rotating case either), but on the set of all approximate ground states. 
These are defined as sequences of (bosonic) A^-particle density matrices 7Ar 
(that is, positive operators on TY^r with trace one) with Tr i^^vTiV ~ NE^^. 
One can then expect that the reduced one-particle density matrix 7^"* of any 
such approximate ground state is a convex combination of GP minimizers, 

7S^^(x,x')« j:A.<Af^(x)<AfP(x')* 

i 

where each (j)^'^ is a GP minimizer, and A, = A^. 

Theorem [2] below states that this is indeed the case. The mathematically 
precise formulation is slightly complicated by the fact that the set of GP 
minimizers is, in general, not countable. 

Let r be the set of all limit points of one-particle density matrices of 
approximate ground states: 

r = I7 : 3 sequence 77V, lim ^Ti HnIn = E'^^{g,n), 

N^oo, Na^g iV 

Since Hq has a compact resolvent by our assumption ([3]) , one easily sees that 
Tr7 = 1 for all 7 E F. Moreover, because of the linearity of the conditions 
in ([Tj), r is clearly convex. 

THEOREM 2. For given value of g > and ft, let T denote the set of all 
limit points of one-particle density matrices of approximate ground states of 
Hjy, defined in ([^. 

(i) r is a compact and convex subset of the set of all trace class opera- 
tors. 

(ii) Let Text C T denote the set of extreme points in P. We have Pext = 

: = E'^^ {g,VL)} , i.e., the extreme points in P are 

given by the rank-one projections onto GP minimizers. 

(iii) For each 7 G P, there is a positive (regular Borel) measure dfi.y, 
supported in Pext; with Jp ^ diJ,^{(p) = 1, such that 

1=1 dfi^{(p) 10) (01 

Text 

where the integral is understood in the weak sense. That is, every 
7 € P is a convex combination of rank-one projections onto GP 
minimizers. 

We remark that item (iii) of Theorem [2] follows from item (ii) by Choquet's 
Theorem |17j . 
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As explained above, Theorem[2]is the natural analogue of ([6]) in the rotat- 
ing case. It can also be interpreted as a rigorous proof of superfluidity. As 
typical for superfluids, angular momentum in rotating systems is acquired in 
terms of quantized vortices. These can be seen by solving the GP equation. 

Theorem [2] also shows the occurrence of spontaneous symmetry break- 
ing. As remarked earlier, axial symmetry of the trap V{x) leads to non- 
uniqueness of the GP minimizer for g large enough [12^ I13j . Uniqueness 
can be restored by perturbing Hq to break the symmetry and favor one of 
the minimizers. This then leads to complete BEC in the usual sense, since 
r contains contains only one element in case the GP functional ([1]) has a 
unique minimizer. 

As in the case of the ground state energy discussed in the previous sub- 
section, the situation is very different for the absolute ground state. The 
set r consists of only one element in this case (namely the minimizer of the 
density matrix functional discussed below, which is unique for any value of 
Q and g). In particular, there is no spontaneous symmetry breaking in the 
absolute ground state. This will be discussed in the next subsection. 

2.3. The Absolute Ground State. Let £'abs(A^5 a, Q) denote the absolute 
ground state energy of H]\f in (j4j), irrespective of symmetry constraints, i.e.. 

Note that necessarily £'abs(A^5 O) ^) ^ Eo{N,a,^}). As is well known, for 
= the two energies are equal. This turns out not to be the case for 
O 7^ 0, in general. 

In the GP limit, the absolute ground state energy, and the corresponding 
one-particle reduced density matrices of approximate ground states, turn 
out to be described by a GP density matrix functional, introduced in |12j . 

£^^[^] = Tr [Ho^] + Ang [ p,{xfd\ . 

Here, 7 is a positive trace class operator on L^(M^), and denotes the den- 
sity of 7, i.e., p-y{x) = 7(x, x). The functional can be shown ^12j to have 
a unique minimizer (under the normalization condition Tr7 = 1), which 
we denote by 7^^. We denote the corresponding energy by E^^{g,Vl) = 

The following Theorem concerning the absolute ground state of Hf^ was 
proved in [13]. 

THEOREM 3. For any fixed g > and Q e M^ 



Ei,^siN,g/N,n) A V ^ (1) DM 



Here, 7^!,^ denotes the one-particle density matrix of any approximate 
(absolute) ground state sequence of Hj^. In other words, the set T defined 
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as in ([7]), but for the absolute ground state, contains only one element, 
namely the unique minimizer of f^^. 

Note that is the restriction of to rank one projections. In the 
case of symmetry breaking (i.e., for g large enough), rank > 2, and 
hence < E'^^ . In particular, in view of Theorems [iHSl the absolute 

and bosonic ground state differ significantly, in general, both in terms of 
their energy and their reduced one-particle density matrix. 

We remark that the results explained in this subsection become physically 
relevant if one considers bosons with internal degrees of freedom. Inter- 
nal degrees of freedom effectively increase the number of allowed symmetry 
classes (see, e.g., jl8j). In particular, if the number of states of the internal 
degrees of freedom of the bosons is greater or equal to the rank of 7^^, 
E^^{g,D.) equals the (bosonic) ground state energy per particle in the GP 
limit. More generally, one can show that in the GP limit the functional £^^, 
when restricted to density matrices of rank at most n, correctly describes 
the ground state energy (and corresponding one-particle density matrix) of 
bosons with n internal states. 

3. Sketch of the Proof of Theorem [T] 

In the following, we shall give a brief outline of the main ideas in the 
proof of Theorem[TJ For details we refer to the original work in [Ij. We shall 
restrict our attention to the appropriate lower bound on the ground state 
energy EQ(N,a,^}). The corresponding upper bound can be obtained via a 
variational argument, as explained in [13|. 

A convenient way to keep track of the bosonic symmetry requirement is 
to work in Fock space. Recall that the bosonic Fock space J- is given by 
T = ©jv>o^^- -'^^ terms of creation and annihilation operators aj and 
on J^, the Hamiltonian can be written as 

^ = Y1 + ^ X] ^i^ki ■ (8) 

j>l ijkl 

Here, we choose the basis in the one-particle space L^(M^) as to diagonalize 
Hq, i.e., Hq = Ylj ^j\'Pj){'Pj\j ^^'^ ^] creates a particle with wavefunction 
ipj, whereas annihilates it. The coefficients Wij^i are given in terms of 
expectation values of Vaix), namely Wijki = {'^i ® ^j\va\^k ® ^i)- 

Note that H in ^ commutes with total particle number operator Ylj>i 
Hence it splits into a direct sum of operators on Ti^ for = 0, 1, . . . . In 
fact, our Hn in (j3D is just the restriction of H to Hn- 

The analysis employed for obtaining a lower bound on the ground state 
energy of H in the sector of N particles consists of two main steps: 

1. Eq. ([8]) is not necessarily well defined. E.g., if Va{x) is the hard-core 
interaction potential (or, more generally, is not integrable), then 
Wijki = 00 for any set of indices. In order to overcome this problem, 
we shall first show that, for a lower bound, one can replace Va{x) by 
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a "soft" and longer ranged potential U (x) (with the same scattering 
length), at the expense of the high-momentum part of the kinetic 
energy. We note that this step is necessary even in the case when 
Va{x) is integrable (and hence ([8]) is well defined) in order to proceed 
with the second step. 
2. After having replaced Va{x) by the softer potential U{x), one then 
shows that it is possible to replace the operators aj and by com- 
plex numbers Zj without changing the ground state energy too much 
[19] . Note that if all the a] and flj in ([8]) are treated as numbers, the 
expression ([8]) looks very similar to the GP energy functional ([I]); in 
fact, it is given by 

(0z|i/o|0z) + ^ / Va{x-y)\Mx)\'^\My)\'^d\d\, 

with (/)z(x) = Y^j Zj(pj{x). 
In the following, we shall explain these two main steps in more detail. 

3.1. Step 1: Generalized Dyson Lemma. The following Lemma can be 
viewed as a generalization of an idea of Dyson [20]. The purpose of the 
lemma is give a lower bound on the interaction potential Va{x) in terms of 
a softer and longer ranged potential U{x), at the expense of some kinetic 
energy (see also [15|). For our purpose, we can only spare the high momen- 
tum part of the kinetic energy, however; the low momentum part is needed 
for the Hq term in the GP functional. 

We thus have to separate the high momentum from the low momentum 
part of the kinetic energy. This can be done in the following way. The proof 
of Lemma [1] is given in [21] . 

Lemma 1. Let Va{x) have scattering length a and range Rq. Let Bji be the 
characteristic function of the ball {x : \x\ < R}. Let < xip) ^ 1; such 
that h{x) = 1 — x{x) is bounded and integrable, 

fnix) = sup \h{x -y) - h{x) \ , 
\y\<R 

and 

wr{x) = -^fnix) / fR{y)d% . 

Then for any e > and any positive radial function U{x) supported in 
Rq 1^ \x\ < R with J U = Air we have the operator inequality 

- yx{p)OR{x)x{p)V + \va{x) > (1 - e)aU {x) - ^wr{x) . (9) 

Here, xip) denotes a multiplication operator in momentum space. Note 
that the operator —V x{p)(^ r{x)x{p)^ can be interpreted as a Laplacian 
that has been localized to the ball of radius R and cut off in momentum 
space. Because of the cut-off, this is not a local operator, however. The 
parameter R is chosen such that a <^ R <^ N^^^^. Note that to leading 
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order in a/ii, the scattering length of 2aU{x) is given in terms of its first 
order Born approximation as (87r)~^a f^s U{x)d?x = a. 

Because of the appearance of the characteristic function 6ji{x) in ([9]), 
Lemma[T]has the fohowing immediate consequence. If yi, . . . , ?/„ are n points 
in whose mutual distance is at least 2R, then 

n n 
i=l i=l ^ 

This bound accomplishes the replacement of the hard interaction potential 
Va{x) by a soft one, at the expense of the high momentum part of the 
kinetic energy. For given configuration of — 1 particles, this estimate 
is applied to the remaining particle. Of course one still has to estimate 
the contribution from configurations where 2 (or more) of the — 1 fixed 
particles are closer together than 2R. This can be achieved by a Feynman- 
Kac integral representation [22] of the ground state. We refer to \T\ for 
details. 

3.2. Step 2: Coherent States. The Fock space J- can be viewed as an 
infinite tensor product of the form T = with spanned by the 

vectors (aj)"|0) for n = 0, 1, . . . . Here, |0) denotes the Fock space vacuum. 

Consider first the case of a single mode, J-i, say. For z € C, a coherent 
state [23] in !Fi is defined by 

\z) = e-l^l'/2+^'^l|0). 

These states span in the whole space J-i. In fact, they satisfy the complete- 
ness relation 

/ dz\z){z\ =1, (10) 

where dz stands for ■7T~^dxdy, and z = x + iy, x, y G M. 

In terms of coherent states, upper and lower symbols of operators can be 
defined. Lower symbols are simply the expectation values of operators in 
coherent states, e.g., (z|a]^|z) = z and (2;|a|a-|^|z) = Upper symbols, 

on the other hand, represent functions of z which, when integrated against 
\z){z\dz over C, yield given operators. For instance, it is not difficult to see 
that = f dzz\z){z\, while a\ai = f dz{\z\'^ — l)\z){z\. Hence, upper and 
lower symbols of are given by z, whereas the lower symbol of aja^^ is |zp 
and the upper symbol is |zp — 1. 

Note that lower symbols yield upper bounds on ground state energies, by 
the variational principle, while upper symbols are useful for lower bounds. 
The difference in the symbols thus quantifies the error one makes in replacing 
the operators a\ and by numbers. In particular, for every quadratic term 
a\a^ a factor —1 has to be taken into account. For this reason, one cannot 
introduce coherent states of all the modes j, but only for a finite number of 
them. 
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In fact, we shall introduce coherent states of all the modes 1 < j < J for 
some J ^ 1. That is, we first write J-" = J-< ® ^>, where J-<^ is spanned by 
the vectors of the form • • • {a\Y-^\^), with Uj G N for 1 < j < J. For 

z = (zi, . . . , zj) G C"^, we introduce the projection operator n(z) on 
given by 

n(z) = \zx ® ■ ■ ■ ® zj){zx ® ■ ■ ■ ® Zj\ . 
Using upper symbols, we can then write the Hamiltonian in ([8]) as 

11= \ dzn(z) (g) /i(z) . 

Here, /i(z) represents the upper symbol of H. Since only the modes 1 < 
j < J have been replaced by numbers, h{z) is an operator on Using the 
completeness property of the coherent states, Eq. (fTOl) . it is then easy to see 
that 

inf spec H > inf inf spec h{z) . 

z 

One then proceeds to show that h{z) ~ S'^^ [(f>z] modulo controllable error 
terms. These error terms are, in fact, operators on ^> which describe both 
the interactions among particles in high modes as well as the interaction 
between particles in modes j < J and j > J. Precise bounds on these terms 
can be found in [Ij. 

4. Sketch of the Proof of Theorem [2] 

In order to obtain information on (approximate) ground states from bounds 
on the energy, one proceeds as follows. One first perturbs the Hamiltonian 
Hn in dl]) by some one-particle perturbation 5, and applies the same per- 
turbation to the GP functional ([T|). One then shows that the result of 
Theorem [1] still holds for the perturbed system. In fact, the proof of Theo- 
rem [1] outlined in the previous section is sufficiently robust in order to easily 
incorporate such a modification. 

Griffiths' argument |24j then implies that, for any 7 G P, and any bounded 
hermitian operator S", 

TV57 > min (0|5|(/>) , (11) 

(f>=(j)G'P 

where the minimum on the right side is taken over all GP minimizers. In- 
equality (jlip is the key to the proof of Theorem [2j The rest follows from 
convexity theory [25], as we shall explain now. 

Recall that an exposed point of a convex set C is an extreme point p with 
the additional property that there is a tangent plane to C containing p but 
no other point of C. Hence, for 7 G P an exposed point, there exists an S 
such that 

Tr 57 < TV 57 for all 7 G P . (12) 
with equality if and only if j = 

It is not very difficult to show that G P for any GP minimizer 

(/>^^. Hence, if we choose 7 in (|12p to be equal to \(f>'^^){(f>'^^\ for the (j)'^^ 
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that minimizes the right side of (lllh for this particular S, the inequahties 
(fTTI) and ([12]) imply that 



and hence there is actually equality in (jl2p . This, in turn, implies that 
J = |(/)*-'^)((^^^|. We have thus shown that all exposed points of T are of 
this form! 

In order to extend this result to all extreme points, now merely exposed 
points, we employ Straszewicz's Theorem L25j, which states that the ex- 
posed points are a dense subset of the extreme points. Strictly speaking, 
this theorem only holds in finite dimensions and not, a priori, in the infinite 
dimensional case under consideration here. However, because of compact- 
ness, the set r is "almost" finite dimensional, and hence the theorem can be 
applied via an approximation argument. We refer again to Ij for details. 



We have presented a rigorous justification of the Gross-Pitaevskii approx- 
imation for sufficiently dilute rotating Bose gases. For large particle number 
N and both Na and of order 1, the ground state of a rotating Bose gas 
is well approximated by the solution to the GP equation. This is true both 
for the energy and the reduced density matrices. In particular, our analysis 
proves the appearance of quantized vortices and the occurrence of sponta- 
neous symmetry breaking in the parameter regime where these phenomena 
can be observed in the GP equation, e.g., for J7 7^ and g large enough. 

We point out that one of the major open problems in this field is the 
validity of the GP equation for rapidly rotating gases, where either |r2| ^ 00 
as ^ cxD (in case the trap potential grows faster than quadratic at infin- 
ity), or Q approaches the trap frequency (for traps that are asymptotically 
quadratic). There is evidence that the GP descriptions breaks down once 
the number of vortices in the system is of the same order as the number 
of particles. Despite recent progress in this direction [26j, a proof of this 
assertion is still lacking. 
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